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We discuss matrix elements of the strangeness changing vector current 
using their relation, due to analyticity, with irK scattering in the P- 
wave. We take into account experimental phase-shift measurements in 
^ ' the elastic channel as well as results, obtained by the LASS collabora- 

tion, on the details of inelastic scattering, which show the dominance of 
two quasi two-body channels at medium energies. The associated form 
^ factors are shown to be completely determined, up to one flavour sym- 

^ ■ metry breaking parameter, imposing boundary conditions at i = from 

chiral and flavour symmetries and at t ^ oo from QCD. We apply the 
Q ' results to the r — > Kt:Vt and r KirnVr amplitudes and compare the 

j-^i former to recent high statistics results from B factories. 
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1 Introduction 

P ' Decays of the r to hadrons with strangeness S = —1 can be used to determine basic 

parameters of the standard model such as Vus and the mass of the strange quark (e.g. [H |2] 
^ I for recent updates). The considerable improvement in statistics brought in by Babar and 

^ ' Belle should translate in the near future into much more precise measurements of matrix 

elements of currents with 5 = — 1 than possible at present. In this paper, we re-consider 
in some detail the relations between the Kit matrix element of the strangeness-changing 
vector current and Ktt scattering in the P-wave. 

We will follow a method applied some time ago by Donoghue, Gasser and Leutwyler 
(DGL [3]) to the vrvr matrix element of the 5 = scalar current (which is not directly 
accessible to experiment since the Higgs boson is not very light). For this method to 
operate, it is necessary that inelastic scattering can be approximated in terms of a finite 
number of two-body or quasi-two-body channels in a sufficiently large energy range E ^ Eq. 
Then, imposing a limited number of constraints at £^ = from chiral symmetry and at 
E » Eq from asymptotic QCD one can determine the form factor from the T-matrix. 
In practice, this is done by solving a set of coupled Muskhelishvili-Omnes (MO) integral 
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equations which are consequences of analyticity properties of the form factors and of time 
reversal invariance. 

Recently, this method was implemented in the case of the K-k scalar form factor [4]: 
in this case, inelasticity in the 5-wave is saturated by the Kr] and K-q' channels. One 
motivation for our interest in the vector form factor is the availability of good experimental 
data from the LASS collaboration on both elastic [5| and inelastic Kti scattering [6l [7l [8] 
in the P-wave. In particular, these works show that in an energy range E ^ 2.5 GeV, 
inelasticity is dominated by two quasi-two-body channels: K*t: and Kp. This makes it 
possible to probe the DGL method by comparing its results with experimental data from 
r decays. The spectral function for r KttVt can be described, of course, by a simple 
Breit-Wigner ansatz in the vicinity of the K*(892) resonance, where scattering is nearly 
perfectly elastic. Away from the resonance, however, this is no longer true. In the energy 
region E « Mk* the spectral function gets dominated by the scalar component [9j (as 
it is less suppressed by phase-space than the vector component). In the energy region 
E » Mk* it is important to correctly treat inelastic effects. Another motivation for 
this work is in view of applications to three body non-leptonic -B-decays B Kmr. In 
kinematical configurations where the kaon and one pion are quasi aligned, factorization 
can presumably be justified |10j and the amplitude gets expressed in terms of Kn vector 
and scalar form factors [Hjl. 

The plan of the paper is as follows. After introducing some notation for the form 
factors involved we discuss the construction of the T-matrix from fits to the experimental 
ttK scattering data. Three-channel unitarity is enforced using the /C-matrix method, and 
proper flavour symmetry and low energy behaviour are enforced starting from a resonance 
chiral Lagrangian. Next, we introduce the MO integral equations satisfied by the form 
factors and discuss their resolution. This necessitates to make a plausible ansatz for the 
T-matrix at high energies. The ansatz determines the number of conditions to be imposed 
in order to determine the form factors from solving the integral equations. We will use 
three conditions at -E = and one asymptotic condition. Finally, the results are displayed 
and compared with recent experimental measurements. 



2 Strangeness changing vector current form factors 

The object which will mainly interest us is the Kir matrix element of the strangeness 
changing vector current 



< K+{pK)\u'y''s\7r^{p^) > 
with t = (pk — Pn)"^ and we will denote the vector form factor as 



H,{t) ^ U{t) = V2f^^-\t) . (2) 



^ In the context of B decays, the usefulness of appealing to descriptions more sophisticated than Breit- 
Wigner combinations for scalar form factors was pointed out in ref. [12] . 
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As we will discuss below Ktt inelastic scattering in the P-wave is dominated by two quasi 
two-body channels: K*tt and pK. This leads us to introduce the associated vector current 
matrix elements 

{p\pv,\)\u^^,s\K-{pk)) = -e^uape*%\)p'^p^j,H^{t) 

(3) 

(we have chosen different signs in the definition of H2 and such that only plus signs 
appear in subsequent equations). In the following, isospin symmetry breaking will be 
neglected. Prom the isospin decomposition of K~^Tr^ it follows that 

(i^+7r0|u7^s|0) = yi([K7r]i|t27^s|0) (4) 

and analogous relations hold for K*tt, Kp. In order to derive the unitarity equations 
satisfied by the form factors Hi^ H2 and -^3 it is convenient to focus on one of the spatial 
components of the vector current, say /U = 3, and go to center-of-mass (CMS) frame of the 
meson pair. In this frame, the current matrix elements introduced above can be expressed 
as follows 

{[K'K]i\u'y^s\0) = V6 cos9oqKAt)Hi{t) (5) 

where qk-k is the modulus of the CMS momentum of the meson pair and Oq is the polar 
angle of the momentum vector with respect to a fixed frame. Similarly, the matrix elements 
which involve one vector meson read 

{[K*{X)TT]i\u-fh\0) = ir]*^ smeoVtqK'nit) H2{t) 

{[p{X)K]i\u-ih\^) = ii]*^ sin 00 Vt qpK{t) Hsit) . (6) 

Here, r/ is an arbitrary phase which can be introduced in the definition of the vector meson 
polarization vector 

e*{X = ±1) 

(^0) 4'o being the polar angles of the vector meson momentum). In the following, we will 
set irj* = 1. We have also taken the following convention for the Levi-Civitta tensor 

£0123 = 1 • (8) 

The explicit dependence on displayed above in eqs. ([5]) ^ indicates that these matrix 
elements concern the angular momentum state J = 1 of the meson pair. 



^ —A cos 6*0 cos (^0 — ^ sin 
— Acos6'osin</)o + icos</)o 
\ Asin^o 



(7) 



3 



3 Kn scattering in P-wave with n-channel unitarity 



3.1 Experimental situation 

Detailed partial-wave analysis of Ktt Kn scattering (in the energy range E < 2.5 GeV) 
have been performed based on high statistics production experiments Kp — > KttN in 
refs. [131 [5j. The LASS collaboration has also performed production experiments K~p — > 
K2-kN and K~p — > K SttN [6l[7l[8] which provide informations on inelastic Kn scattering. 
These, however, are not as precise as for elastic scattering and mainly concern resonance 
properties in the various partial- waves. In the 5-wave, Kn scattering is elastic to a good 
approximation up to the Krj' threshold|13l [5]. indicating that K-q' is the main inelastic 
channel. Inelasticity in the P-wave starts at somewhat smaller energy than in the S- 
wave and the results of refs. [6l suggest that K*n is the main inelastic channel, as it is 
enhanced by two resonances i^*(1410) and K*{1680) (following the nomenclature of the 
PDG jUj ). The latter resonance was observed to couple essentially to three channels: 
Kn, K*n as well as Kp [14]. This indicates that Kp should be also taken into account as a 
significant inelastic channel in the P-wave although its coupling to K*(1410) was found to 
be small. The experimental results on the branching ratios of the K*(1410) and K*(1680) 
are collected in table [H 

The amplitude Kn — > Kt] was studied in ref.[8] and found to be very small in the P- 
wave and not to display any resonant effect. The coupling of Kr]' to resonances with J^"-^ = 
1 is proportional to the sine of the mixing angle and thus should also be suppressed. 
These results suggest that a plausible description of Kn scattering in the P-wave with n- 
channel unitarity can be performed (up to energies E < 2.5 GeV) by retaining as inelastic 
channels the two quasi two-body channels with one vector meson: K*n and Kp. 





Kn K*n Kp 


K*(1410) 
K*(1680) 


6.6 ± 1 > 40 < 7 
38.7 ± 2.5 29.9l|^ 31.4l|j' 



Table 1: Experimental decay branching ratios (in percent units) for the resonances 
i^*(1410) and K*(1680) quoted in the PDG. 



3.2 Lagrangian for vector resonances and pseudo-scalar mesons 

It is useful to express the resonance contribution to scattering in terms of coupling constants 
which are known in other contexts and it is also useful to impose chiral constraints at low 
energy. For this reason, let us start from the following Lagrangian which was used in the 
context of the chiral expansion[T5l [T6]. It includes the nonet of the light vector mesons 
encoded in a matrix and the light pseudoscalars encoded in a chiral matrix and 
involves two coupling constants gy and ay, 

£«=£(J)±£«±£« (9) 
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with 



= ]^av e^^^P-ti {V^{u,, Vl) (10) 

and V^i, = V ^Vy — V yV^. In ref.[16], for example, the following values are quoted for the 
coupling constants (based on the extended NJL model) 

gv ^ 0.083, av ^ 0.25 (11) 

which should serve as a guide as to the orders of magnitudes. Let us denote an excited 

(n) 

vector resonance by , we can write down the following coupling terms 

cP = ^gvin)tr{V^:^^[u„Uy]) 

4") = iav(n) e'^'^'^^tr (Vi^Huy, V}) . (12) 

These terms do not involve the quark mass matrix and, therefore, have exact 5^7(3) flavour 
symmetry. Experimental results on K* resonances indicate that flavour symmetry can 
sometimes be substantially broken, as in the case of the ii'*(1410) (see tabled]). We do 
not try to write down all possible Lagrangian terms which break flavour symmetry but 
eventually will implement such effects at the level of the fits. We also do not consider 
explicitly the possibility of many more terms which involve further derivatives acting on 
the vector or on the chiral fields. Again, such terms which give rise to polynomial energy 
dependence, may be implemented phenomenologically as "background" contributions in 
the fits. 



3.3 Resonance contributions to Kn scattering 

Let us recall the constraints introduced by the conservation of parity. The Kit system, at 
first, in the P-wave has parity —1. The action of the parity operator on a state formed 
from a vector meson and a pseudo-scalar meson involves the helicity[T7]. 

V\V{X)P>j={-y-'\Vi-X)P>j . (13) 

For J = 1, the following combination is the only one which has negative parity 

^- = ^ {\Vil)P >j=i -\Vi-l)P >j=i) . (14) 

We will need the partial-wave expansion of the scattering amplitudes in a situation where 
the initial state CMS momentum has polar angles 4'o and the final state has polar 
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angles 9, 4> [^7] 

(y (A)P|r|i^7r) = levr ^ (2 J + 1) {qvp qx^Y 



J,M 

{V{X)P\T-'\K7T)Vll^{(t), e)Vi^,icl>o, Oo) 

{V{X)P\T\V'{X')P) = IdTT ^ (2 J + 1) {qvP qvp] 

J,M 

{ViX)P\T-'\V\X')P)Vli^{^,6)ViryicPo,eo) . 



(15) 



We have factored out explicitly the angular momentum barrier factors in the definition of 
the partial-wave T-matrix elements: this is necessary for non-diagonal matrix elements in 
order to ensure good analytical behaviour for the partial- wave T-matrix [18| . For simplicity, 
we also define the diagonal elements in the same manner e.g. 



{Ktt\T\K7t) = levr ^ (2 J + 1) (qx^) 



2J 



X 



J,M 

{K7T\T-'\K7T)VllQ{<f>,e)Vi,^o{(t>o,eo) . (16) 

Let us number the three relevant channels in our construction as 

1 — > Ktt, 2 — > K*TT, 3 — > Kp, (17) 

Using the Lagrangian introduced above (fTOl) . (fT2]l it is not difficult to compute the res- 
onance contributions to the various T-matrix elements. The resulting amplitude can be 
written in a compact form which displays the usual resonance structure , 

rpi _ V- gr{n,i)grin,j) 

^ij,r - 2^ ^^^> 

with 




9r{n,3) = -— . (19) 

Vlovr r-n- 

The relations between channel 2 and channel 3 matrix elements are consequences of exact 
fiavour symmetry as implemented in the Lagrangians (fTOl) . (fT2]) . The energy dependence 
of these effective coupling constants is reliable only for small values of s ( s << M^, ) 
where chiral symmetry constraints (as encoded in the construction of the Lagrangian) are 
relevant. In practice, we will implement a simple cutoff function by replacing 



(l + 6, 



l + bn 



(20) 



AS 



in eqs. (fT9|l . This cutoff function has the correct behaviour at small s and effectively 
replaces -y/s by the resonance mass M„ otherwise. 
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3.4 ii'-matrix fits to the experimental data 



We may implement n-channel unitarity in a simple way, using the -fC-matrix method 
(e.g. |19|). Starting from a real, symmetric X-matrix we define the J = 1 T-matrix in 
the following way (suppressing the J superscript), 

r= (l + iKSQ2)-ii^ (21) 
where and E are diagonal matrices 

^Q^\i3 — ^ijOij '^ij — ^ij~/= (22) 

Qi being the CMS momentum in channel i. Indeed, it not difficult to verify that the 
5-matrix , which is defined as follows, 

S =l + 2iV^QTQV^ (23) 

is unitary 

SS^ = 1 (24) 
and encodes the proper J = 1 angular momentum barrier factors. 



3.4.1 A first simple fit 



Mn gv{n) (7v{n) 

0.894 0.0714 0.0989 

1.719 0.0103 -0.1857 

2.247 0.0126 -0.3069 



(25) 



The most detailed experimental results concern the Ktt elastic channel. We used the results 
of ref. [Ij on the amplitudes K~7r~^ — > K~7r~^ in the energy range E < 2.5 GeV and those 
of ref. [13] on the isospin 3/2 amplitude. Let us first perform a simple fit including just 
three resonances in the X-matrix, i.e., we using the formulas (fTSl) (I19p . We 

find a qualitatively acceptable result (see figs.[T]) with the following resonance parameters. 



(masses are in GeV and the coupling constants are dimensionless) . In this fit we have taken 
6„ = oo (see ([20]) ) i.e. ^/s is replaced by M„. This fit does find a resonance corresponding 
to /C*(1680) in addition to the K*{SS^2) but prefers to locate the third resonance at a 
higher energy rather than at 1.4 GeV. Obviously, though, the energy region between 1 and 
1.5 GeV is where the data is not very well described (the overall /dof = 7.4). 



3.4.2 More sophisticated fits 

Including just one more resonance does not improve the situation very much. In order 
to obtain significantly better fits we will include a fourth resonance in addition to non- 
resonant background terms. There are many physical sources for such terms. We have seen. 
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3 resonances fit 







K-matrix — 
Aston ' — 
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Simplest 3 resonance fit 
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E.K (GeV) 



Figure 1: K-matrix fit to the Kir Kn data of ref./^ with three resonances, ap is the 
moduhis of the K~Tr^ i^~7r+ amphtude and $p is the phase. 



for instance, that the couphng constants may carry energy dependence. In addition, we 
expect contributions associated with the left-hand cut, i.e. arising from meson exchanges 
in the crossed channels. We parametrize such contributions in the following, simplistic, 
way 

j^back _ "^l* 



K 



back _ g(Q2 + ass) 

12 - l + ,3 



^13 =YT7^- (26) 

For the other -fC-matrix elements we do not introduce any background dependence. We 
have tried many possibilities but, clearly, the amount of experimental data is insufficient 
for probing in detail all the matrix elements of the K matrix, so one must make admittedly 
arbitrary simplifying assumptions. We found that setting non-zero background terms for 
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^11 Ki2, is the most economic way (in terms of number of parameters) to achieve a 
good fit. In addition to the elastic Kir data, we try to reproduce the constraints on the 
inelastic channels in the resonance regions (see tabled]). These data imply strong flavour 
symmetry violation in the region of the K*(14:W) resonance. We account for this fact by 
relaxing the symmetry relation gr{n,3) = —gr{n,2) (see eqs. (fT9]) ) for n = 2. Instead of 
this relation, we suppress the coupling to channel 3 (i.e. Kp ) by simply setting 

5r(n,3)=0 (27) 

when n = 2. 

Altogether, this fit contains 16 parameters and allows for a rather satisfactory fit to 
the Aston et al. data, which has a x^/dof ~ 1.8. We note that the fourth resonance 
present in this fit effectively acts essentially as an additional source of background at lower 
energies: it should not necessarily be interpreted as a true physical resonance (in fact 
no corresponding resonance is listed in the PDG). The numerical results for the best fit 
parameters are collected in table [2l 



n 


Mn 


gv{n) 


o-y(n) 




1 


0.8962 


0.72820 10" 


-i 


0.26080 


3.3906 


2 


1.3789 


0.5252310-^ 


^2 


-0.56075 


18.373 


3 


1.7300 


0.6936510-^ 


^2 


-0.21774 


-9.2048 


4 


2.2739 


0.1204410-^ 


-2 


-0.29360 


2.8318 



Table 2: Results for the fit parameters concerning the four resonances and the background 
as described in sec. \3.4.2\ 

One observes that the values of gv{^) and cry^l) are in reasonable agreement with the 
ENJL predictions from ref. [16]. Fig. [2] shows the comparison of the fit results with the 
experimental data. The improvement with respect to the result of the simpler fit shown in 
fig. [U is obvious, particularly in the energy region [1,1.4] GeV. Moreover this fit is able to 
reproduce, qualitatively at least, the experimental results concerning the inelastic channels 
recalled in table [H This is illustrated in fig. [3] which shows the moduli of the transition 
5-matrix elements and, for the elastic channel 1 — 5ii (the cross-sections are proportional 
to the squares of these quantities). Indeed, the figure shows that in the energy region 
corresponding to the K*{1AW) the matrix element 5i2 shows a clear resonance peak while 
5i3 shows no peak. In the energy region of the K*(1680) resonance these two matrix 
elements display peaks of the same height (this was imposed as a constraint in the fit) and 
the peak in the elastic channel is slightly higher. 

In the very small s region finally, our T-matrix is expected to be qualitatively reasonable 
but certainly not very accurate because of the lack of constraints from experimental data 
in the threshold region. The ttK scattering length, for instance, is found to be 

a^^ = 0.025 (28) 

which has the correct order of magnitude but is slightly larger than the result from ChPT 
(at order /)t20j: aj^^ = 0.018. 
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4 resonances+background fit 







K-matrix 
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0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



E.K (GeV) 

Figure 2: K-matrix fit to the Kn Kn data of ref.|^ with four resonances plus back- 
ground, as described in sec. \3.4.2\ 



4 The Ktt vector form factor from n-channel MO equations 
4.1 Unitarity equations 

The Muskhelishvili-Omnes equations derive from the fact that, firstly, each of the form 
factors Hi{t), H2{t) and {{^{t) satisfy an unsubtracted dispersion relation because they 
are analytic functions in the variable t except for a right-hand cut E| [H] and they decrease 
for t — > oo faster than 1/t. Secondly, one can express the imaginary parts of the form 
factors in terms of T-matrix elements using time-reversal invariance. Let us briefly repeat 
the derivation in the case of a scalar operator O which is T-invariant 

TOT-^ = O . (29) 

^In reality, the form factors H2, H3 could display anomalous thresholds since they involve an unstable 
particle. We ignore this possibility in our analysis. 
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4 resonances+background fit 




We consider the imaginary part of the matrix element between the vacuum and a state 
\m > which we take to be an "out" state 

2ilm < mlOlO >=< mlOlO > - < mlOlO >* (30) 

out out out 

Using the fact that T is antiunitary and transforms an "out" sate into an "in" state, it is 
not difficult to transform the second term on the right-hand side in (f3Q]l and recover the 
usual unitarity expression for the imaginary part 

Im < m\O\0 >= ^ E < , > ■ (31) 

out 2 — out 

n 

In reality, we will use operators which carry space-time indices and thus are not T-invariant. 
It is however easy to see that the T variation of the operator is compensated by the T- 
variation of the momenta pi in the state \m > such that eq. ([3T]) remains valid for the 
form factors. For a given energy of the state \m >out a finite number of states |n >out 
contribute to the sum in the right-hand side (Watson's theorem [22] follows if the energy 
is lower than the first inelastic threshold). If we truncate the summation in eq. (fSTI) and 
insert the imaginary parts into the dispersion relations we obtain a closed set of MO integral 
equations for the form factors. 



4.2 Application to Kn, K*7t, pK vector form factors 

We start from the general expression (fSTl) . apply it to the operator n7^s and retain three 
states in the sum: Kvr, K*7r, pK . We then use the formulas ^ ([6]) and the partial- 
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wave expansions (fTSll of the T-matrix elements and compute the phase-space integrals. 
As expected, only the J = 1 partial-wave contributes and, after a small calculation the 
unitarity equations for Hi{t), H2{t) and H-^{t) are obtained. They can be written as 
follows, in matrix form 




H2 



(32) 



where T is the 3x3 J = 1 T-matrix, the diagonal matrices and S are defined in eq. ([22 
and r is also a diagonal matrix 



diag (1, Vi) . 



(33) 



The MO integral equation set then derives by combining eqs. ([32]l with the unsubtracted 
dispersion relations satisfied by the form factors 



/ Hi{t) 



1 




(34) 



We recall here that the following matrix plays a role when discussing existence and multi- 
plicity of the solutions to the MO equations[23l [2^ 

(35) 



S = 1 + liT-^TQ'^T.T 



This matrix differs from the 5-matrix , as defined in eq. (|23]) , but the determinants of the 
two matrices are equal 

det 5 = det 5 . (36) 



4.3 Asymptotic conditions on the T-matrix 

The MO equations obeyed by the form factors are coupled, homogeneous singular integral 
equations with a kernel linear in the T-matrix. The mathematical properties of such equa- 
tions are exposed in Muskhelishvili's book[23]. In particular, the number of independent 
solutions N is given by the index of the integral operator which can be expressed in terms 
of the sum of the eigenphases 5i{t) of the iS-matrix. For an n x n S'-matrix, 

j2Mc^)-5,m=N^ ■ (37) 
1 

In order to determine the form factors from the integral equations one must therefore 
impose independent conditions. In our case one has n = 3 and the ^-matrix has been 
constrained from experimental input up to Eq ~ 2.5 GeV. We make a key assumption here 
that Eq is sufficiently large that the asymptotic regime for the T-matrix sets up for values 
of E not much larger than Eq. At E = Eq the sum of the eigenphases is J2 ^ii^o) — 3.57r. 
We therefore expect the index N to be either 3 or 4. We choose to adopt an asymptotic 
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ui-Q')^. ^ "r (38) 



condition which has an index = 4. This will enable us to impose four conditions on the 
form factor. Three of these use the values at the origin of the form factors, Hi{0). We 
will discuss below how these values are constrained by experiment. As a fourth condition, 
we can enforce the behaviour of at infinity. In QCD, ignoring flavour symmetry 

breaking, one should have|25l [26] 

We do not attempt to reproduce the logarithmic running of Us, and actually use eg. ((38]) 
with a constant value Us = 0.2. The condition ((381) is then easily implemented in the form 
of a sum rule 

1 rco 

- / dt'lmHiit') = 16^^/2a,F2 . (39) 

In the region t > to one must use a parametrization which respects the unitarity of the 
5-matrix . For this purpose, we write the 5-matrix in exponential form, 

S = exp{2iH) (40) 

where H is a real, symmetric matrix. We then define the interpolation on the matrix H 
such that 

lim Hu{t), Hi3{t), H2z{t) = 0, 

r— >oo 

lim Hii{t) + H22{t) + H33{t) = Att , (41) 

I— »oo 

using the simple interpolating function 

H,,{^)(l-(^^yy . (42) 

The parameters Oij and Pij are determined such as to ensure continuity of Hij(t) and its 
first derivative at t = to. Consistent with the assumption made above on the setting of 
the asymptotic regime, the parameter 7 must be larger than one: in practice, we will take 
7 = 2. The condition on the trace of H still leaves some freedom as to the behaviour 
of each diagonal element. We will consider a plausible scenario in which the asymptotic 
eigenphases satisfy 

5i(oo)=27r, 62(00) = 53(00) = TT (43) 
where 61 is the largest eigenphase at t = to- The three eigenphases are shown in fig. [H 



4.4 Conditions at t = 

In this section we discuss the experimental constraints on the form factor components at 
t = 0. The component Hi(0) = /+(0), firstly, is well known from ChPT to very close to 1, 
the result at O(p^) was computed in ref. |27j, 

/+(0) = 0.977 (ChPT 0(/)) . (44) 
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S-matrix eigenphases 
2tt I 1 1 1 




Vi (GeV) 



Figure 4: Eigenphases of the S-matrix: the left part of the figure shows the region deter- 
mined from experiment, the right part (colored in yellow) shows the asymptotic interpola- 
tion region. The asymptotic condition here is J2 ^iioo) = 47r. 



We note that ChPT computations at order have been performed [28] but the coupling 
constants involved are not yet well known. The remaining two form factor components 
involve vector mesons. We will argue that their values at t = can still be evaluated 
reasonably well by appealing, in addition to chiral expansion arguments, to a leading large 
Nc approximation. In the chiral limit, flavour symmetry is exact, and the following relation 
holds between a charged current matrix element and an electromagnetic current one, 

(K* + |u7^.|vr°) = ^(/,+ |j^^^,|vr+) . (45) 

This relation implies that one can relate i?2(0) and -^3(0) to the radiative width of the 
charged p meson, 

1^,(0,1 Ha(0)|=(fl^)\ (46, 

Using the experimental value [l^ of the radiative decay width 

rp+^^^+ = 68 ± 7 keV (47) 

gives 

\H2{0)\ = \H-i{0)\ = 1.54 ± 0.08 GeV"^ . (48) 

This method does not fix the sign of -^2(0). For this purpose, we can appeal to a simple 
vector-dominance picture applied, e.g., to the vertex function {'y\u^f^d\7r~) . In such a 
picture -^2(0) gets related to the ABJ anomaly [29], 

ft(0) = -if3(0) = ^^.1.50 GeV-'. (49) 
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We observe that the absolute value of -^2(0) in the VMD model is in rather good agreement 
with the one deduced from experiment. 

We can try to refine these estimates by taking into account the breaking of flavour 
symmetry to first order in the quark masses. For this purpose, let us write down an 
effective Lagrangian, 

= w{ Vtr(F'^K,/f}) 

+ ^tr({y^n'^}{x+,/f}) 

+ ^tr([y^x+]K/f])}. (50) 



We have used the same notation as in ref. [16] for the first term in this Lagrangian. Flavour 
symmetry breaking effects can be encoded into three independent Lagrangian terms: this 
holds in the leading large N^. approximation (since multiple trace terms are suppressed). 
We can then constrain the coupling constants hy, a,b and c by considering the radiative 
decays p"*" — > 7r+7, K*^ — > ^"'"7, and K*^ K^^. The corresponding amplitudes, 
computed from the Lagrangian (|50p . have the expressions 



- ^+7) = /^y^^ h'Jf = hy + ''^h (51) 



and 



AiK*-^ K^,) = _ W2 - 3c) 

^ " 3Fk ^ 3F^ ^ ' 

A{K*' ^ i^°7) = =^h^yff - . (52) 

^ " 3Fk ^ 3F^ ^ ' 



The following experimental results are available [14] for the K* mesons radiative decays, 

T{K*+ ivr+7) = 116 ± 12 keV, 

r(K* ^ K^j) = 50 ± 5 keV , (53) 

which, together with the result ([47l) allows one to determine three coupling-constants, 

hlf^ = 0.0356 ± 0.0018, 
5 = 0.0010 ±0.0015, 

c = 0.0032 ± 0.0013 . (54) 
Let us now compute our form factor components from the Lagrangian ([SOll . obtaining 
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These expressions show that i?2(0) and -^3(0) depend on the symmetry breaking parameter 
a which is left undetermined by the analysis of radiative vector meson decays. Numerically, 
using the values ([54ll we find 



We expect the parameter a to have the same order of magnitude as the other two symmetry 
breaking parameters, i.e. \a\ < 10~^. 

5 Solutions and comparisons with experimental results 

5.1 Solving the MO equations 

The integral equations ((32]) ([341) can be solved numerically by discretizing the integrals. 
This must be done in a way which guarantees a precise evaluation of the principal-value 
integrals. For this purpose, we use expansions over Legendre polynomials and exact ex- 
pressions for the principal- value integrals of these. More details can be found in ref. |24| . 
One then arrives at a set of M x M homogeneous linear ordinary equations. The Muskhe- 
lishvili conditions on the number of solutions imply that the determinant of the system 
must vanish (if it does not, then the only solution would be the trivial identically vanishing 
one) and the associated matrix must have N (with = 4 in our case) zero eigenvalues. 
In practice, because of discretization and roundoff errors no eigenvalue vanishes exactly, 
but one does have eigenvalues which are very small in magnitude. A numerically stable 
way to proceed is to enlarge the system to an (M + N) x M one by adding constraints 
on the form factors as equations and then solve the new system using the singular value 
decomposition method. We obtained precise results with M > 100 and we used values of 
M up to 400. Several correctness and accuracy tests can be performed. In particular, while 
exact solutions are not known for general T-matrices, the value of the determinant of the 
matrix formed from independent solutions can be expressed in analytical form[23l [24] 
in terms of 5'-matrix elements. 

5.2 Results for r decays 

The vector form factor f+{t) can be probed using the r decay mode r KnUr. The 
energy distribution of the Kn pair has the following expression which involves f+{t) as 
well as the scalar form factor fo{t) 



H2{0) = 1.41 ± 0.09 - 65.4a GeV~^ 
HsiO) = -1.34 ± 0.07 - 65.4 a GeV" 



-1 



(56) 




(57) 



( 




2^ \ 
J t 



\U{t)\' + 



3(m2 



K 




\fo{t)\' 
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where the definition of /o, in terms of the form factors introduced in eq. ^ is 



Mt) = V2 



(5^ 



The contribution of /o is kinematically suppressed except for very small values of t. Let us 
begin by fixing the symmetry breaking parameter a from the integrated width. Assuming 
exact isospin symmetry one has 

3 

Txtt = 3rj^-^o = . (59) 



The following value for the rate is quoted by the PDG|14| 

rP^G ^ ^^2,.f> ± 0.05) 10"^ . (60) 

The most recent results from the Babar and Belle collaborations[30l [31]tend, however, to 
point towards a slightly smaller value 

^Babar ^ ^^2A% ± 0.009 ± 0.054) 10-^ 

^Beiie ^ ^^2.12 ± 0.006 ± 0.039) 10"^ (61) 

(assuming isospin symmetry). These results can be reproduced with our vector form factoiEl 
by having the parameter a in the range 

a = (-7.01°:^) 10-3 . (62) 

The central value corresponds to the PDG result and the value a = —9.0 10"^ reproduces 
Belle's central figure for the rate. The energy distribution in the decay r — > K-ki^t has been 
measured for the first time by the Aleph collaboration[33|- The Belle collaboration has 
recently measured the distribution in energy for the decay r~ — > KgiT'Vr with considerably 
higher statistics [3l] (approximately a thousand times larger). The Babar collaboration 
has also presented results for the energy distribution for t~ K' [30] but the data 

corrected for background are not publicly available. Data from Babar on r~ — > Kgir'UT 
have been analyzed in a thesis [34] but the results have not yet been published. 

The result from our construction is compared with the result of Bell^ in fig. [5] for 
two values of a. In addition to the (dominant) contribution from the vector form factor, 
we have also included the contribution from the scalar form factor which we computed 
following ref. The parameter a affects the size of the K*{892) peak (which controls 
essentially the value of the integrated decay rate) and also the shape of the form factor in 
the inelastic region. We find that it is possible to reproduce both the integrated decay rate 
and the shape of the energy distribution above the K* mass reasonably well. The solutions 
in the region \/t > 1.4 GeV are sensitive to the assumptions made on the 5-matrix in the 



■'ill applications to r decays we use for Vus the value such that Vusf+{0) = 0.2167 quoted by the 
averaging group [32J- 

^ The data shown have been corrected for background but not for acceptance, which we have assumed 
to be approximately energy independent. 
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Vi (GeV) 

Figure 5: Results for the energy dependence of the r decay rate r KnUr using the vector 
form factor as discussed in the text for two values of the symmetry breaking parameter a 
compared with the result from Belle ISDj . The contribution from the scalar form factor fo 
is shown dash-dot curve. 

asymptotic domain. For instance, if we choose an asymptotic condition with = 3 rather 
than = 4 we cannot get agreement with experiment when \/i > lA GeV by varying a. 

It is instructive to compare also our result for the vector form factor with that obtained 
by Belle from a fit to their data. This comparison avoids the problem of the acceptance 
(which is taken into account in their fit) but one must keep in mind that there is some 
model dependence in this determination (a model independent separation of the vector 
and scalar form factors requires to analyze angular distributions) . In ref . [31] five different 
fits, using descriptions of the form factors /+, /o a la Kiihn and Santamaria [3^, have 
been performed. In fig. [6] we compare our result for the modulus of the vector form factor 
with the experimental determination based on the fit which includes two vector resonances: 
K*{892), K*(1410) for the vector form factor and one scalar resonance: Kg (800) for the 
scalar form factor. For clarity, we have presented our model's result as one curve corre- 
sponding to the best fit to the LASS data but, obviously, one should keep in mind the 
uncertainties of these data, which are sometimes sizable, e.g., in the 1.3 GeV region (see 

fig. El). 

This comparison deserves a comment. In the K*{892) region one expects to observe 
some isospin symmetry breaking effects since the Kn system in r ^ Kt:Vt- is in a charged 
state while the experimental data from LASS concern Kvr in a neutral state. Fig. [6] 
shows a visible difference concerning the width which is narrower for the K*^ than for the 
Somewhat surprisingly, no difference is seen concerning the mass. In fact, the mass 
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t2 




Figure 6: The yellow band represents the modulus of the vector form factor determined 
from a four parameter fit in ref. j3Jf It is compared with our determination for two values 
of the symmetry breaking parameter a. 



resulting from Belle's fit: M^.+ = 895.47 ±0.20 ±0.44 ±0.59 differs by about 4 MeV from 
the mass quoted in the PDG: Mj^*+ = 891.66±0.26 which is based on hadronic production 
experiments. This shift in the K*~^ mass was reported earlier by CLEO [36] but remains 
to be confirmed. It could be a similar effect to the one observed for the p meson mass. 

The model also generates predictions for r decays into K*7r and Kp via vector current. 
The energy distribution of the decay width for K*tt reads, in terms of the form factor H2 

(63) 

and an exactly similar expression holds for Kp in terms of H^. The results from our 
construction are plotted in fig. El The figure shows that the K*(1410) appears very clearly 
in the K*it channel. The r decay into Kp (via vector current) is comparatively strongly 
suppressed. These features reflect the constraints which we have used in the construction 
of the T-matrix based on the LASS experiments. The results for the integrated rates are 

R{t K*TTVr)v = (1.37 ± 0.02) 10"^ 

R{t Kpvr)v = (4.5 ± 3.0) 10"^ (64) 
where the errors are estimated by simply varying the parameter a. Our result for R{t 



dV 



K* 



327r3 
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K*iTi^T-)v is seen to agree with the one quoted by the Aleph collaboration [33j, 

{lA+li tH) X 10-3 . (65) 

which, however, is not very precise. 




1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 



Vt (GeV) 

Figure 7: Results for the energy distribution of the r decay width into K*tt and Kp 
(rescaled by a factor of 10) via vector current from our model, for two values of the 
symmetry breaking parameter a. 



5.3 Some results at t = 

For values of t near t = 0, experimental results on can be obtained from Ki^ decays 

and several new experiments have been performed recently. One usually defines the slope 
and curvature parameters from the Taylor expansion 

hit) = /+(0) f 1 + X'^ + \ A"^ + ..\ . (66) 

The results which we get for A' and A" are 

A' = (26.05iJ]|i) 10-3, A" = (1.291^4) 10-' (67) 

where the errors, again, simply reflect the range of variation of the parameter a (the 
lower values correspond to a = — 910~3). They compare reasonably well with the average 
(performed in ref. j32j ) over the recent experiments 

A' = (24.8 ± 1.1) 10-3, A'' = (1.61 ± 0.45) 10-^ . (68) 
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The results ([67)1 can also be compared with the predictions recently presented in ref. [9]n . 

A' = 26.2 10~^ A" = 1.37 10^^ ( ref. [9]) . (69) 



5.4 Some remarks on the asymptotic region 

When i ^ oo, the MO equations driven by the T-matrix behaving as discussed in sec. 14.31 
can be shown to imply that tf+{t) goes to a constant. We have constrained this constant, 
via eq. (f39l) . in order to correctly reproduce QCD in an average sense. Fig. [8] displays the 
real part of the product obtained from the numerical solution of the equations. The 

figure also shows the asymptotic QCD expectation (analytically continued to the timelike 
region). One sees that t/+(t) indeed goes to a constant, but the asymptotic behaviour sets 
in at fairly large values of t. A similar feature was observed in the case of the pion vector 
form factor which was discussed in ref. [37| . 




Vt (GeV) 

Figure 8: Real part of the product t/+(t) resulting from our solution of the MO equations 
compared with asymptotic QCD prediction. 

Another remark concerns the phase of f+{t) (let us denote it as (p-^^{t)). It is often 
assumed that (/>+(*) should go to vr as t goes to infinity (e.g. [38]), while compatibility with 
asymptotic QCD simply implies that it should go to vr modulo 27r. Let us recall that if the 
phase of /+(t) goes to vr at infinity, /+ can be expressed in terms of its phase as a minimal 
Omnes representation. In our construction, we find that the behaviour of the phase at 
infinity depends on the value of a. This is illustrated in fig. [H For the central value of a 

^We converted the numerical values taking into account that the expansion formula l[66|l is used in 
ref. [9] with = 138 MeV rather than = 139.57 MeV. 
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the phase actually goes to Svr at infinity. When a gets slightly smaller, a transition occurs 
and the phase goes to vr. Again here, the asymptotic behaviour is reached for rather large 
values of t. 



a 



a 



a 



7.010-3 
7.510-3 
9.010-3 



3^ - 



27r - 



TT - 








1 



10 



Vt {GeV) 



Figure 9: Illustration of the changing behaviour of the phase of f+{t) at large t when the 
parameter a varies. 



6 Conclusions 

In this work, we have considered matrix elements between light states of the stangeness 
changing vector current. Such quantities are now becoming accessible to experiment, with 
rather good precision in a rather large energy range from hadronic r decays. There are 
possible applications to certain three-body B decays as well. From a theoretical point 
of view, one can establish relations between two-body form factors and the scattering 
matrix not only at very low energy in the elastic scattering region (Watson's theorem) 
but also at somewhat higher energies provided that inelastic scattering is dominated by 
two-body channels. This was shown to be the case for ttK scattering in the P-wave by the 
LASS collaboration[6l [3, [5] in the energy range E < 2.5 GeV. Similar properties hold for 
TTTT scattering in the 5- wave and also for ttK scattering in the 5- wave and were used to 
construct scalar form factors^, Those results, however, have not been compared with 
experimental data. 

We have constructed a three channel T-matrix, which satisfies unitarity, from fits to 
the large amount of data on ttK ttK scattering in the P-wave and satisfying the exper- 
imental constraints on inelastic scattering. In the region E > 2.5 GeV, one must make a 
plausible guess. We impose a smooth interpolation such that the MO operator has index 
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= 4. The three form factors can then be determined by solving the system of MO 
equations and applying four constraints. We have used the three values at the origin, 
which we have argued to be determined, making use of chiral and flavour symmetries, up 
to one symmetry breaking parameter a. As a fourth constraint, we used asymptotic QCD 
behaviour for the form factor /_|_. 

The result of this construction for the form factor f+{t) compares reasonably well with 
the recent determination by the Belle collaboration. The agreement in the inelastic region 
is not quite trivial to achieve. Indeed, the -fC*(1680) resonance appears as a very large 
effect in ttK ttK scattering while it is suppressed in the form factor. Varying the 
parameter a outside of the range allowed from the integrated r — s- Kttz/^ rate also destroys 
the agreement in the energy distribution. The result in the domain \/t > 1.4 GeV is 
also sensitive to how the T-matrix behaves above 2.5 GeV. In the region of the K*{892) 
resonance, there are some discrepancies. While these could be due to isospin breaking 
effects, they seem to concern the width rather than the mass of the K*(892), which is 
unexpected. Finally, we made predictions for the total rates and the energy distributions 
for the vector current contribution to the r decays r — > K*t:Vt- and r — > pKur. 
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